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Chapter 1 
Introduction 


The electrical, macroscopic performance of isotropic material can generally 
be described through their constitutive scalar parameters, permittivity and 
permeability which are symbolically represented by e and //, respectively. 
These parameter* relate the electric and magnetic flux densities to the 
electric and magnetic fields through the following relationships: 


D = eE 


( 1 . 1 ) 


and 


B = fiH. (1.2) 

It is through these parameters that the interaction of electromagnetic waves 
with material can be quantized in terms of reflection and transmission co- 
efficients, and propagation and attenuation factors. 

These parameters can be complex which indicates the loss component 
of the material. The loss is represented in the imaginary part of the param- 
eters as e = s' — je" = e 0 e r and /z = - j/z" = /z 0 Mr where the subscripts o 

and r refer to the free space and relative values. When the loss mechanism 
is conductive in nature, e" is replaced with cr/o; where cr is the conductivity 
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of the material and u) is the radian frequency. Often, the loss is stated as 
loss tangents being tan£ e = e M / 1' and tan£ m = 

The acquisition of constitutive material parameters are often obtained 
from reflection and transmission measurements using a coaxial or wave- 
guide fixture with an internal material sample. The parameters are then 
recovered from exact expressions relating the measured quantities to the 
parameters. Obviously, the number of unknowns determine the number of 
independent measurements required. The most convenient measurements 
are reflection and transmission for such closed (nonradiating) fixtures. Ac- 
curate parameters are difficult to obtain when the material can not be 
shaped into the desired shapes resulting in undesirable air gaps between 
the sample and fixture. 

Alternate techniques have been examined to acquire the necessary num- 
ber of independent measurements. The goal has been to develop a technique 
to conveniently measure material samples at high temperatures (greater 
than 800° F) over a broadband of frequencies (2-18 GHz) to determine the 
constitutive material parameters. The simplest approach examined was to 
measure the reflected field as a function of frequency from a waveguide aper- 
ture terminated with a material coated planar conductor. This approach 
would be very desirable since planar material samples can readily be fab- 
ricated. Although the measured reflected field varied between samples and 
sample temperature, there were deficiencies. Since only one measurement 
was possible (reflection), another independent quantity had to be gener- 
ated. One approach was to differentiate or integrate around the desired 
frequency of interest to form another entity. This approach had reasonable 
success [1] where the measurement entailed a free space quantity and a good 
analytical representation of the scattered fields was known. This approach 
was not successful in this case primarily because a sufficient model of the 
reflected field was not readily available. Another potential limitation for 
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this approach is that the second entity may be too similar to the measured 
quantity to form a second independent entity. 

Other alternate techniques entailed the swept frequency measurement 
of scattered fields in a parallel plate or ground plane fixture as described 
in [2]. Comparable scattering measurements for small dielectric spheres 
in free space to obtain the constitutive parameters have had good success 
[3,4]. The parallel plate or ground plane fixture offers a means to heat a 
sample to the desired temperature and perform broadband measurements 
simultaneously. A major drawback to either of these two techniques would 
be the fabrication of cylindrical or hemispherical material samples. In this 
report, parameter determination using the parallel plate and ground plane 
fixtures are discussed with results presented. 


3 



Chapter 2 

Scattering Measurement 
Parameter Extraction 


The method of parameter extraction from measurements is dependent upon 
the fixture used. Conventional coaxial and rectangular waveguide fixture 
uses reflection and transmission measurements in deterministic expressions 
to arrive at the desired parameters [5]. These expressions are available since 
the field behavior can be described with one mode. Performing broadband 
measurements involves many modes in general and a more complex pro- 
cedure to extract the desired parameters. One approach involved the free 
space measurement of a dielectric sphere [4]. For electrically small spheres 
(a < .2A) only two modes (dipole moments) are required to adequately de- 
scribe the scattered field. As the sphere becomes electrically larger, more 
modes are required and accordingly more independent measurements (one 
measurement for each mode). 

An alternate approach was taken to require only two independent mea- 
surements to arrive at the desired parameters. A search algorithm was em- 
ployed to search for the unknown constitutive parameters using the mea- 
sured data and analytical representation. The search algorithm adapted 
here was the Newton-Raphson technique. The major disadvantage for this 
approach is the sensitivity of the algorithm to converge due to noise in the 
data. 
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The Newton-Raphson technique iteratively solves for the unknown pa- 
rameters by forcing two convenient functionals to zero. The functionals, 
Fo i , for each frequency are defined to be 

F 0 - Mq — Co = 0 (2-1) 

and 

F l =M 1 -C 1 = 0 (2.2) 

where the M 0i i and Co.i are the independent measured and calculated re- 
sponses for the sample geometry. The calculated responses are dependent 
upon the parameters e and fj, which are iteratively varied until some mini- 
mum error criterion has been achieved. The unknown parameters are iter- 
ated with the following process from some initial guess 

£t+i = £{ + A £,• (2-3) 


and 


/*;+ i — + A (Xi 


where 


Ae 



and 


A /i — 


IT ill 
en 




( 2 . 4 ) 


(2.5) 


( 2 . 6 ) 


and 
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BFi BF\ 
Be Bp 


(2.7) 


Reliable results can be obtained if accurate measurements are made. 
The collection of these measurements are facilitated with the HP-8510B 
Network Analyzer. This computer controlled device can readily perform 
reflection and transmission measurements between 2 and 18 GHz. A very 
desirable feature of swept frequency measurements, which the HP-8510B 
incorporates through software, is the time gating of undesired signals from a 
desired signal. This time gating eliminates error terms that originate from 
other scattering centers located away from the desired scattering center. 
Such potential error signals are due to scattering from fixture terminations 
which are not included in the analytical model. 


The measurements obtained have to be calibrated before they can be 
used in the search algorithm. The calibration entails the minimization 
of scattering which is not sample dependent and the frequency correction 
of the system response. The measurements required are the frequency re- 
sponses for when the unknown sample is in the fixture and when it is absent 
from the fixture along with the frequency response for a known sample (ref- 
erence) which is usually a perfectly conducting sample. This measurement 
data is then used in the following formula to obtain calibrated data 

CSD = ( mrd-mbd ) crd (2 ' 8) 

where CSD is calibrated sample data, MSD is measured sample data, MBD 
is measured background data, MRD measured reference data and CRD is 
calculated (exact) reference data. The subtraction process in the above 
equation minimizes the influence of other scattering centers that are not 
sample dependent. The system response is normalized from the measure- 
ment through the scaling of the sample response with the measured and 
calculated responses for the reference. 
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Chapter 3 

Measurement Fixtures 


Of the possible fixtures that would allow the samples to be heated without 
disturbing the measurements and having an exact analytical scattered field 
representation for the sample are the parallel plate waveguide and shaped 
ground plane fixtures. Such fixtures were fabricated based upon electro- 
magnetic constraints given in [2] and room temperature measurements were 
performed to evaluate these fixtures. 

The consideration of sample size for these fixtures is that the sample 
can be heated uniformly and the electrical thickness is sufficient to permit 
a significant field interaction to occur. If the sample is supported on a 
conductive substrate and the coating thickness is too thin, the scattered 
field will be influenced more by the substrate parameter (its radius) and 
not by the material’s constitutive parameters. 

I Parallel Plate Fixture 

The parallel plate fixture is simply a parallel plate waveguide which will 
allow an electric field to propagate between the plates to illuminate a sam- 
ple. Figure 3.1 illustrates the fixture which was built using a wood frame 
and aluminium sheeting to form the conductive sheets. The fixture has 
two tapered sections attached to an interior region which the sample is 
measured in. The tapered sections allow the electric field to be directed to 
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and from the interior region for monostatic and bistatic reflection measure- 
ments. The vertices of the tapered feed sections have a coaxial port which 
permits a connection to the HP-8510B Network Analyzer. It is extremely 
important to physically position the coaxial port as close to the vertex as 
possible to minimize an undesired reflection coming from the vertex. The 
spacing between the parallel plates was chosen to allow only the dominant 
mode to propagate in this structure between 2 and 18 GHz. 

The canonical sample shape is a circular cylinder which should maintain 
a tight physical contact between the ends of the cylinder and the aluminium 
plates. Any gap between the plates and sample will result in error in the 
measurement. For even sample heating, the material under test could be 
coated onto a hollow tube with internal electric heating coils. 



Figure 3.1: Parallel plate test fixture. 

The exact scattered fields for the dielectric circular cylinder or a coated 
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perfectly conducting circular cylinder are readily calculated by eigenfunc- 
tion expressions. The geometry for a coated perfectly conducting circular 
cylinder is shown in Figure 3.2. The scattered far field for plane wave 
illumination is given by [6] 

/ 2 _ P ~i k P 00 

(3-D 


where a n are ratios of cylindrical Bessel functions, J n , N n and H ^ and 
given by 


a n = 


c J' n {kb) - ZMM) 
£n ffn^'(kb) - Z„ffi !> (kb) 


(3.2) 


and 


" Vr Uk,<,)N„(hb) - N„(ha)J„(k t b) ' ' 

where t} T = is the relative wave impedance of the material coating and 
fci = ^/e r fi T k and k is the free space wave number. The parameter e n is 
equal to unity for n = 0 and 2 for n ^ 0. 

These above expressions reduce to simpler forms for special cases. The 
quantity Y n is equal to 

v ’ JJMb) 

for a homogeneous cylinder. Additionally, for the perfectly conducting 
cylinder, Y n = oc. 

The 2-18 GHz swept frequency response for a 2” diameter perfectly 
conducting cylinder is shown in Figure 3.3 where the solid curve is the 
monostatic (reflection) case and the dashed is the 90° bistatic reflection 
case. These two cases are the nominal choices to acquire two independent 
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Figure 3.2: Geometry for a coated cylinder. 

measurements. Note that these responses are ideal for calibration refer- 
ences since they are fairly constant. Figure 3.4 illustrates the backscatter 
frequency response for a homogeneous 2” dielectric cylinder (solid line) and 
a 2” diameter perfectly conducting cylinder coated with a .07” thick layer 
of dielectric. In both cases, the constitutive parameters are t T = 24— j2 and 
Hr = 2 — j. 2. Illustrated in this figure is the effect of loss on the scattered 
field for a given amount of material. For the homogeneous cylinder, the 
scattered field becomes solely dependent upon the reflection from the illu- 
minated face of the cylinder as the frequency is increased. This behaviour 
is not the same for the hemisphere case to follow next. In this geometry, a 
creeping wave will exist along with the specular field. 

II Ground Plane Fixture 

The ground plane fixture shown in Figure 3.5 simulates an infinite con- 
ducting surface for which the scattered fields from a hemisphere can be 
calculated. The ogival shape controls the scattered field from the edges of 
the fixture. The dominant scattered fields from the fixture would be the 
tips which are low in value compared to those scattered from a straight 
edge at normal incidence. This fixture allows several different types of 
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Figure 3.3: Swept frequency response for a 2” diameter perfectly conducting 
circular cylinder. Solid line - backscatter, dashed line - 90° bistatic. 

independent measurements to be performed. The measurement depends 
upon the polarization of the incident field, the elevation angle of the inci- 
dent field above the ground plane and the bistatic angle between the source 
and receiver. The most advantageous pair of measurements are when the 
source and receiver are located at the same point (backscatter) and both 
polarizations are measured for the same elevation angle. The source fields 
for this fixture were generated using a broadband 2-18 GHz AEL [7] horn 
connected to the HP-8510B Network Analyzer. 

The canonical sample shape is a hemisphere which should maintain a 
tight physical contact to the ground plane fixture. Any gap between the 
ground plane and sample will result in error in the measurement. For 
even sample heating, the material under test could be coated onto a hollow 
hemisphere with internal electric heating coils or a gas burner. 

The geometry of a coated dielectric hemisphere over an infinite ground 
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Figure 3.4: Swept frequency backscatter response for circular cylinders with 
e T = 24 — j 2 and p r = 2 — j. 2. Solid line - 2” diameter homogeneous 
cylinder, dashed line - 2” diameter perfectly conducting cylinder with .07” 
thick dielectric coating. 

plane is shown in Figure 3.6. The exact scattered far field for plane wave 
illumination can be calculated by using image theory which results in the 
sum of the backscattered and bistatic scattered fields for a sphere in free 
space. The scattered field for a hemisphere on a ground plane can be written 
in the following form for the two principal polarizations 


Ee = E'g vhtTt {0 = 0) + E^ € (8 = 20 o ) (3.4) 

and 

E* = E ,ph ' re {9 - 0 ) - E; phere (6 = 20 o ) ( 3 . 5 ) 


where Eg p £ ere is the scattered field for a sphere in free space at a bistatic 
angle of 0 O . 
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Figure 3.5: Ogival shaped ground plane fixture. 


The scattered field expressions for a sphere in free space are given by 


[3] 


E‘ e phere {0 ) = - 


T? OO 

E 


and 


kr 


E 0 


n~ 1 


b n sin 0P R (cos 0) — c n 


P n{ COS0) 


6in 0 


*?*■"(») = E 


n— 1 


A ^n( CO ) • nplf/ 

b n . c n sin vP R (cos u j 

sin v 


where P^costf) and P 1/ (cos^) are Legendre polynomials and 

BnjA^ - r, r7 X n J*) 


b n = b° 


A n {B n K-r, r2 X n TZ) 


and 


c n — c. 


( X J* - VrtAJ*) 


(3.6) 


(3.7) 


(3.8) 


(3.9) 
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where 77,2 = \[^ is the relative wave impedance of the material coating. 
The symbols f>°, c°, A n , B n ,X n and F^’ 2,3,4 are ratios of spherical Bessel 
functions j n and h n given by 


2n + 1 
n(n 4- 1) 


[kbj n {kb)]' 

kbh { n\kb)] n 


2 n + 1 jn(kb) 

n(* + l)hg\kb)' 


[JcbjJJeb)l 

kbj n (kb) 

B _ [kbhW(kb)] 
kbhL 2 \kb) 

x _ Min(M)]' 

" Mjn(M) 


F2 = Fi - 


Tl = Fi - 


F 3 = Fi - 


*?rl 

[kiaj n (k! 

o)T_ 

[fc 2 ahj, 2 l(fc 2 a) 

1r2 

[k 1 aj n (k 1 

«)] 

fc 2 ahi, 2 \fc 2 a) 

In 

[k\aj n (ki 

o)T_ 

Jfc 2 ah(, 2 )(fc 2 a)j 

lr2 

[*i ajjjkx 

«)] | 

fc 2 ahl 2 i(fc 2 a)] 

Vrl 

[k 1 aj n (k 1 

«)] 

k2ah^\k 2 a) 

Vr2 

[kiaj n (ki 

«)]' 

|fc 2 a/in 2 ^(fc 2 a) 


-F 4 


F 2 

x n 1 


and 



7n \ha-j n {k x a)} 

k 2 ah^\k 2 a)j 

lr 2 [k\aj n {k y a)}' 

|fc 2 ohl 2 *(fc 2 a)j 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 
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where ij rl = is the relative wave impedance and ki = y/F^jl^k in the 

homogeneous sphere region, and k 2 = k with k being the free space 

wave number. The functions F^’ 2,3,4 are ratios of spherical Bessel functions 
given by 


F_ x 


= 1 . - 



[fc 1 oj n (fe iq))' h^jjeib) 

M^b) [jfe ia fcW(fcia)] M 


[Ajajn^a)]' 


[Mjn(M))' | 

kiah!^\kia) 


F 3 _ 1 _ JnjMa) fe l &/t n 2) ( fe l fe ) 

[kibj n (kib)]' h£\kia) 


and 


j n (fc 1 q) hWfab) 

jn(kib) h^ik.a)' 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


These above expressions also reduce to simpler forms for special cases. 
The FjJ’ 2,3 ' 4 quantities reduce to F^’ 2 ’ 3,4 for a perfectly conducting coated 
sphere, and F x — F 3 and F 3 = F 4 for a homogeneous dielectic sphere. 
Additionally, for the perfectly conducting sphere, b n = 6° and c n = c°. 


The 2-18 GHz swept frequency backscatter response for a 2” diameter 
perfectly conducting hemisphere with an incidence angle of 30° from graz- 
ing is shown in Figure 3.7 where the solid curve is the theta polarized return 
case and the dashed is the phi polarized return case. These two cases are 
the nominal choices to acquire two independent measurements. Note that 
the responses at an elevation angle of 30° are reasonable for calibration 
references since they are fairly smooth for both polarizations. Figure 3.8 
illustrates the backscatter, theta polarized frequency response for a homo- 
geneous 2” dielectric hemisphere (solid line) and a 2” diameter perfectly 
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Figure 3.6: Geometry for a coated hemisphere. 

conducting hemisphere coated with a .07” thick layer of dielectric. In both 
cases, the constitutive parameters are t r = 24 — j 2 and fi r = 2 — y.2. Sim- 
ilar comments can be made in regard to loss and material thickness for 
the hemispherical geometries as were made to the circular cylinder results 
shown in Figure 3.4. Also note that the calculations shown in Figure 3.8 
were stopped at approximately 15 GHz due to the recursively generated 
Bessel functions used in the eigenfunction solution. Above this frequency 
for the chosen parameters, the Bessel function subroutines failed. For fur- 
ther calculations, the Bessel function values would have to be generated 
asymptotically. 

Ill Sample Heating 

Of the two fixtures, the ogival ground plane fixture would be the simplest 
to uniformily heat the material sample to the required temperatures. This 
conclusion is formed since the parallel plate fixture has more metal to heat 
to the desired temperature and keep heated to overcome the effective heat 
sink losses. Practical measurements with the ogival plate fixture would 
require that the heating be accomplished from the shadowed side of the 
ground plane and a metal curtain be formed around the perimeter of the 
ground plane to shield the heat source. The parallel plate fixture inherently 
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Figure 3.7: Swept frequency response for a 2” diameter perfectly conducting 
hemisphere with the incident field 30° from grazing incidence. Solid line - 
theta polarizaton, dashed line - phi polarization. 


has a high degree of shielding from potential external disturbances. In both 
fixture designs, it is required that the samples be in good physical contact 
with the metal sides of the fixtures. 
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Figure 3.8: Swept frequency, theta polarized backscatter response for hemi- 
spheres with t T — 24 — j 2 and fJ, T = 2 — j. 2. Solid line - 2” diameter 
homogeneous hemisphere, dashed line - 2” diameter perfectly conducting 
hemisphere with .07” thick dielectric coating. 
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Chapter 4 


Example Measurements and 
Parameter Extraction 


The examples that follow have been chosen to illustrate the concepts in- 
volved. The sample geometries measured were perfectly conducting and 
low loss dielectrics. Both fixtures were used to obtain the necessary mea- 
surements to extract the constitutive parameters for a lossy dielectric. All 
scattering calculations are normalized to echo width for the parallel plate 
measurements and echo area for the ground plane measurements. Echo 
width and echo area are respectively defined as 


aw = 2 it lim p 

p—+oo 


E i 


and 


a a — 4rr lim r 2 

r— * OO 


E‘ 2 

Te* 


(4.1) 


(4.2) 


The calibration technique of the measured data is very important. Equa- 
tion 2.8 involves subtraction and multiplicative processes. The subtraction 
is to remove erroneous scattering terms present in the fixture. Figure 4.1 
illustrates the backscatter terms present for the parallel plate waveguide in 
terms of a time signature. The waveform was generated by transforming 
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the swept frequency response using a fast Fourier transform when a 2” di- 
ameter disk was in the fixture. Figure 4.2 is a similar signature formed by 
subtracting the measurement with the disk absent in the fixture from the 
measurement with the disk present in the fixture. Note that many terms 
disappear but some do not since the fixture was slightly disturbed inserting 
and removing the sample. To remove the other error terms, soft time gating 
has to be employed by either the HP-8510B Network Analyzer or during 
the calibration process. It was found that better results were obtained by 
oversampling (smaller delta frequency increments) in the frequency domain 
and performing a software time gate on the data in the calibration process. 
A better signal to noise ratio can be obtained in the time domain by over 
sampling since the noise is spread throughout a larger time window. 


- 25 . - 20 . - 15 . - 10 . -5 0 5 . 10 . 15 . 20 . 25 . 



Figure 4.1: Scattering time signature of parallel plate fixture illustrating 
the scattering terms present for a reflection measurement. 
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Figure 4.2: Scattering time signature of parallel plate fixture illustrating 
the scattering terms present for a reflection measurement after background 
subtraction. 

The performance of the parallel plate fixture was evaluated by measur- 
ing a known sample. Figure 4.3 compares the measured and calculated re- 
sponses for a 3” diameter, perfectly conducting, .313” thick disk at bistatic 
angles of 0° (backscatter) and 110°. Figure 4.4 shows the same measurement 
conditions for a 2” diameter Teflon disk. Figure 4.5 shows the extracted 
constitutive parameters as a function of frequency. The maximum variation 
around the mean value is approximately 2 %. 

The performance of the ground plane fixture was evaluated by measur- 
ing a known sample. Figure 4.6 compares the measured and calculated 
responses for a 2” diameter, perfectly conducting hemisphere at an eleva- 
tion angle of 30° from the ground plane for both principle polarizations. 
The performance of the system is good between 4 and 11 GHz. Figure 
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4.7 shows the same measurement for a 2” diameter Plexiglas hemisphere. 
Figure 4.8 shows the extracted constitutive parameters for this Plexiglas 
hemisphere as a function of frequency. As with the parallel plate fixture, 
the maximum variation around the mean value is approximately 2 %. 

Other measurements were performed to observe the influence of a small 
gap between the ground plane and test samples. This gap would be desir- 
able to concentrate the heat on the sample by providing an insulating air 
barrier to the metal ground plane. The measurements indicated substantial 
variations from the ideal situation for both incident polarizations such that 
no constitutive parameter extraction would be possible using the existing 
analytical model. 
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Figure 4.3: Parallel plate swept frequency results for a 3” diameter perfectly 
conducting disk: Solid - measured, dashed - calculated, a) Backscatter. b) 
bistatic, angle of 110°. 
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Figure 4.4: Measured parallel plate 6wept frequency results for a 2” diam- 
eter Teflon disk. Solid - backscatter, dashed - bistatic, angle of 110°. 
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Figure 4.5: Extracted constitutive parameters for a 2” diameter Teflon disk. 
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Figure 4.6: Ground plane swept frequency results for a 2” diameter per- 
fectly conducting hemisphere, 30° from grazing: Solid - measured, dashed 
- calculated, a) Theta polarization, b) Phi polarization. 
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Figure 4.7: Measured ground plane swept frequency results for a 2” di- 
ameter Plexiglas hemisphere, 30° from grazing. Solid -theta polarization, 
dashed - phi polarization. 
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Figure 4.8: Extracted constitutive parameters for a 2” diameter Plexigl 
hemisphere. 





Chapter 5 

Extraction Technique 
Sensitivity 


The ability to extract the correct constitutive parameters from material 
measurements depends upon how well the fixture and sample geometry 
simulate the analytical model. Some differences between the physical and 
theorectical models can be tolerated. Errors originating from scattering due 
to fixture discontinuities can be minimized through processing measured 
data if the discontinuities occurred in a different physical location from 
where the sample was located. Errors originating from sample geometry 
differences such as incorrect dimensions can not be corrected. 

An informative study to acquire tolerance limits on the extracted con- 
stitutive parameters is to extract these parameters from calculated swept 
frequency data assuming one of the geometry parameters is slightly differ- 
ent. Figure 5.1 illustrates the extracted parameters, normalized to the free 
space values, for a homogeneous dielectric, 2” diameter hemisphere with 
material parameters of e r = 2.6 — jO. and fx r = 1. — y0.. The incident field 
was theta polarized, 30° from grazing. The curves shown are for an assumed 
diameter of 1.96” and 2.04” for the perfectly conducting hemisphere. The 
maximum variation of the scattered field for these two cases was approxi- 
mately 5 dB. The nominal measurement variation has been approximately 
3 dB. 
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The influence of measurement processing is also of interest. Figure 5.2 
shows a result of different smoothing of the backscattered swept frequency 
field for a 2” diameter Teflon disk in the parallel plate fixture as shown in 
Figure 4.4. The extracted relative permittivity is slightly dependent upon 
the degree of smoothing. Smoothing data with a moving average reduces 
noise and other scattering error terms. Note that a moving average is 
actually a discrete convolution and in the time domain it represents itself 
as a time gating window to remove terms outside the time region of interest. 

Finally, the ability to extract large relative constitutive parameters 
through a search algorithm becomes more difficult as the the value be- 
comes larger. The required accuracy of the measurements becomes very 
severe. Also as the loss becomes large, the functional dependence of the 
scattered field simplifies as shown in Figure 3.4 for the lossy homogeneous 
circular cylinder. This situation limits the extraction of parameters since 
all the scattered fields are basically described by a reflected field which is 
dependent upon only the ratio of e and fi. 
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Figure 5.1: Normalized extracted constitutive parameters from calculated 
scattered field values. Assumed diameter for the perfectly conducting hemi- 
sphere: Solid line - 2 % smaller, dashed line - 2 % larger. 
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Figure 5.2: Backscattered results for a 2” diameter Teflon disk in the paral- 
lel plate fixture. Normalized extracted permittivity: Solid line - 31 points, 
dashed line - 71 points. 
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Chapter 6 

An Alternate Fixture 


The fixtures built and tested for parameter extraction were designed for 
broadband measurements using proven scattering techniques with a single 
fixture. Fabricating material samples in the required geometries demands 
care but such geometry sample shapes have been made in the past including 
spherical shapes. To achieve more accurate measurements than with these 
“free” space scattering approaches, the measurements have to be taken us- 
ing a “non-radiating” fixture. The most common fixture incorporates a 
waveguide. However, the requirement of generating broadband measure- 
ments with a single fixture is then relaxed. 

The conventional technique of waveguide reflection and transmission 
measurements can provide good signal to noise measurements. The stan- 
dard sample approach is to have a sample which would completely fill the 
rectangular waveguide The requirement of the sample completely filling 
the guide is based upon the equations used to extract the constitutive pa- 
rameter values from the measurements. It is proposed to eliminate the 
requirement of having the sample completely fill the waveguide interior as 
shown in Figure 6.1. The sample would only partially fill the waveguide 
interior in one dimension. The proper set of equations would have to be 
developed for this geometry and a comparable Newton-Raphson technique 
could be used to extract the constitutive parameters. The sample geometry 
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has the convenience of being in a more readily fabricated shape and may be 
desirable since it is a planar coating which is closer to the final application 
form of the material. Such a fixture also has the advantage of providing 
high temperature measurements by placing a portion of the waveguide fix- 
ture in an oven. The actual design for a high temperature fixture would be 
different than the one shown in Figure 6.1 due to metal fatigue, oxidation 
and expansion considerations. Presently the sample holder concept would 
have to be changed since potential error terms can result from scattering 
due to the joints between the sample holder and the remaining fixture. 
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Figure 6.1: An alternate fixture using a rectangular waveguide fixture and 
a planar fixture sample. 
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Chapter 7 
Conclusions 


The wideband constitutive parameters were obtained for two generic geome- 
tries from scattered field measurements. A 2 % parameter variation was ob- 
served in the extracted parameters. The ability to acquire accurate enough 
measurements is the major electrical limitation for the presented technique. 
Software processing enhances the accuracy of the measurements by reducing 
the effect of fixture related error terms. An advantage of this approach over 
traditional waveguide techniques for material parameter measurements is 
that the material samples are free to expand for a given temperature. The 
actual high temperature sample dimensions could then be estimated from 
the room temperature dimensions using known coefficients of expansion. 
The ogival shaped ground plane fixture is more desirable over the parallel 
plate fixture for high temperature considerations of providing uniform heat 
to a small area with smaller fixture error terms to calibrate out. 

To achieve more accurate measurements, a closed, nonradiating fixture 
may be considered as shown in Figure 6.1. Heating the fixture and sample 
becomes simpler. The sample geometry is also much easier to fabricate with 
a reduction in sample tolerances as compared to the traditional waveguide 
approach of filling the complete cross section with material. 
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